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ABSTRACT 

Context. Arakelian 564 (Ark 564, z = 0.0247) is an X-ray-bright narrow-line Seyfert 1 galaxy. By using advanced X-ray timing 
techniques, an excess of “delayed” emission in the hard X-ray band (4 - 7.5 keV) following about 1000 seconds after “flaring” light 
in the soft X-ray band (0.4 - 1 keV) was recently detected. 

Aims. We report on the X-ray spectral analysis of eight XMM-Newton and one Suzaku observation of Ark 564. Our aim is to 
characterise the X-ray spectral properties of the source in the light of these recently reported results. 

Methods. High-resolution spectroscopy was performed with the RGS in the soft X-ray band, while broad-band spectroscopy was 
performed with the EPIC-pn and XIS/PIN instruments. We analysed time-averaged, flux-selected, and time-resolved spectra. 

Results. Despite the strong variability in flux during our observational campaign, the broad-band spectral shape of Ark 564 does not 
vary dramatically and can be reproduced either by a superposition of a power law and a blackbody emission or by a Comptonized 
power-law emission model. High-resolution spectroscopy revealed ionised gas along the line of sight at the systemic redshift of the 
source, with a low column density (N/j ~ 10^‘ cm“^) and a range of ionisation states (-0.8 < log(^/erg cm s“') < 2.4). Broad-band 
spectroscopy revealed a very steep intrinsic continuum (photon index T ~ 2.6) and a rather weak emission feature in the iron K band 
(EW ~ 150 eV); modelling this feature with a reflection component requires highly ionised gas, log(^/erg cm s“*) > 3.5. A reflection- 
dominated or an absorption-dominated model are similarly able to well reproduce the time-averaged data from a statistical point of 
view, in both cases requiring contrived geometries and/or unlikely physical parameters. Einally, through time-resolved analysis we 
spectroscopically identified the “delayed” emission as a spectral hardening above ~ 4 keV; the most likely interpretation for this 
component is a reprocessing of the “flaring” light by gas located at 10-100 rg from the central supermassive black hole that is so hot 
that it can Compton-upscatter the flaring intrinsic continuum emission. 

Key words. Galaxies: active - Galaxies: Seyfert - Galaxies: individuals: Ark 564 - X-rays: individuals: Ark 564 - Black hole 
physics 


1. Introduction 


Arakelian 564 (Ark 564) is one of the X-ray-brightest AGN 
of the sky, with a 2-10 keV flux / 2 _io ~ 10“'^ erg s“' 
cm“^. It is classified as na rrow-line Sey fert 1 (NLSl), with a 
FWHM(H/3)~750 km s ' (IStiroel Il99(ll) . Like many NLSls, 
Ark 564 displays an e xtremely variable X-ray flux and a 
very steep spectral shape (iBrandt et alJll994l : lTurner et al.ll2001l : 
lEdelson et alJ200^ . Possible complexities in its X-ray spectrum 


related to the reprocessing of the primar y continuum emissio n 
onto ionised m atter h av e been reported b y lVaughan et al. ( 199^, 


[Turner et al.l d 19991) . IComastri et alJ (l200llk Vignali et akl 


(l2004l). Matsumoto et alJ ( 20041), [Panadakis et alJ (l2007l) . 


iDewangan et al.l (12007 ). Smith et S (12008 ). 


iLegg et ^ (1201 2l) exploited a long XMM-Newton observa¬ 
tional campaign to perform a detailed timing analysis of the X- 
ray properties of Ark 564. The main result was the discovery of 
an excess of emission in the hard (4-7.5 keV) band following 
soft (0.4-1 keV) flares that appeared ~ 1000 s after the flares 
and lasted ~ 1500 s. This is consistent with light reverberation 
on some scattering material located ~ lOOr^ from the continuum 
source (rg = GMbhIc^). 
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The main goal of this work is characterising the X-ray spec- 
tral properties of Ark 564 in the light of the results reported by 
iLegg et ^ (l2012h . The paper is organised as follows; the obser¬ 
vations are presented in Sect. |2] the spectral analysis is reported 
in Sect. [3] the results are summarised in Sect.|4l and conclusions 
are drawn in Sect. |5] 


2. Observations and data reduction 

In the following subsections we report on the data reduction of 
eight XMM-Newton and one Suzaku X-ray observations of Ark 
564. 

2.1. XMM-Newton observations 

Ark 564 was observed by XMM-Newton eight times between 
May and July 2011. Data were processed using HEAsoft v.6. 12 
and SAS v. 11.0.0. 

The reflection grating spectrometer (RGS) data were pro¬ 
cessed using calibration files generated in June 2012. Back¬ 
ground light curves were generated and inspected for the pres¬ 
ence of flares: a uniform threshold of 0.1 ct s“' for the back¬ 
ground rate was then used to generate good time interval ta¬ 
bles for each observation. Net exposure times after background 
flaring removal are reported in Table [T] First-order spectra and 
response matrices were extracted using the rgsproc task us¬ 
ing source coordinates as given by NED (RA=340.663939, 
dec=-i-29.725364) as input to prevent possible offsets between 
different observations. While there is considerable flux variabil¬ 
ity between the different observations, the spectral shape is found 
not to vary, therefore we proceeded with stacking the eight obser¬ 
vations using the task rgscombine. We obtained two high signal- 
to-noise spectra with ~ 412 ks of exposure each, and a total of 
~ 10^ counts for the combination of the two instruments, RGSl 
and RGS2. Source spectra were rebinned by a factor of 3. 

All the EPIC pn observations were performed in small win¬ 
dow mode and used the thin optical blocking filter. The raw 
ODE pn event files were filtered, and we retained only best- 
quality (FLAG=0) single and double (PATTERN < 4) events. 
Source events were extracted from circular regions with a 36" ra¬ 
dius, while background events were extracted from two boxes of 
61" X 100" and 51" x51" in size. Background light curves were 
then extracted for single events alone in the 10-12 keV band 
for each observation. After inspecting these light curves, we re¬ 
moved strong flaring background time intervals using a conser¬ 
vative threshold of 0.1 ct s *. The net exposure times after the 
flaring background removal is reported in Table [1] These values 
also account for the small window mode live time of ~ 71%. 

The SAS task epatplot was then used to estimate the frac¬ 
tion of pile-up affecting the data. By using the standard energy 
interval 0.5 - 2 keV to estimate the pile-up fraction, excellent 
numerical results were obtained, with the single or double pat¬ 
tern observed distributions deviating by less than 1% from the 
predicted ones. However, a visual inspection of the pattern dis¬ 
tribution revealed significant deviations both at energies below 
1 keV (with an observed excess of single events) and at ener¬ 
gies above 1 keV (with an observed excess of double events), 
the two effects compensating in the numerical estimate for the 
pile-up. We then re-estimated the pile-up fraction by focusing 
on the 1-9 keV band. Using this band as a reference, we found 
non-negligible deviations (of the order of 3 -5%) of the observed 
pattern distribution from the predicted one in half of the obser¬ 
vations. The cores of the PSF were then excised using a radius 


proportional to the amount of pile-up in the data, until we were 
able to recover the predicted pattern distribution within a 1 - 2% 
uncertainty over the whole energy range used in this work, that 
is, 0.4 - 10 keV. Excision radii are reported in the last column of 
Tabled Redistribution matrices (RMF) and ancillary response 
files (ARE) were generated for each observation at the position 
of the source using tasks rmfgen and arfgen , taking into account 
the excision of the PSF core. 

The EPIC-MOS data were heavily affected by pile-up and so 
were not useable. 

2.2. Suzaku observation 

Suzaku (iMitsuda et al.ll200^ carries four co-aligned telescopes 
that focus X-rays onto a suite of CCD cameras comp rising the X- 
ray Imaging Spectrometer (XIS iKovama et all2007[) . XIS CCDs 
0, 2, 3 are front-illuminated (FI) and XIS 1 is back-illuminated. 
XIS 1 has an enhanced soft-band response relative to the FI 
units, but has a lower effective area at 6 keV and a higher back¬ 
ground level at high energies. XIS 2 failed on 2006 Novem¬ 
ber 9 and hence was not used in our analysis. Suzaku also 
carries a non-imaging , collimated Hard X-ray Detector (HXD, 
iTakahashi et al.l2007l) whose PIN detector provides useful AGN 
data typically over the range 15-70keV. 

A Suzaku observation of Ark 564 was made in 2007 June 
26-28 (OBSID 702117010), and the data were reduced using 
v6.16 of HEAsoft. We screened the XIS and PIN events to 
exclude data taken during passage through the South Atlantic 
Anomaly, starting and ending within 500 seconds of entry or 
exit. In addition, we excluded data with an Earth elevation angle 
smaller than 5°. Ark 564 was observed at the nominal centre po¬ 
sition for the HXD. FI CCDs used 3x3 and 5x5 edit-modes, with 
normal clocking mode. Good events (with grades 0, 2, 3, 4, and 
6) were selected and hot or flickering pixels removed using the 
SISCLEAN script. The spaced-row charge injection (SCI) was 
used. The net exposure times were 70 ks per XIS unit and 81 ks 
for the PIN. XIS products were extracted from circular regions 
of 2.9' radius while background spectra were extracted from a 
region of the same size, offset from the source and avoiding the 
calibration sources in the chip corners. 

For PIN analysis we used the model “D” background (Q). 
We checked the correctness of the model adopted by compar¬ 
ing it to the Earth-occulted data, finding an overestimation of 
the NXB of the level of 11.5%; this factor was taken into ac¬ 
count into the subsequent analysis. Following standard practice, 
the time filter resulting from screening the observational data 
was applied to the background events model. The ftool hxd- 
piNXBPi was used to create a PIN background spectrum from the 
screened background data, t aking the shape and flux level of the 
cosmic X-ray background (lBoldti[T987t iGruber et al.l|1999ll in 
the Suzaku PIN field of view into account. We used the response 
file ae_hxd_pinhxnome3_20080129.rsp for spectral fitting. 

During the Suzaku observation Ark 564 was found to have 
source count rates 2.132 + 0.006 (summed XIS 0, 3 over 0.75- 
lOkeV) and 0.010 +0.002 (PIN over 15-50keV) ct/s. The back¬ 
ground was 1% of the total XIS count rate in the full band. The 
source comprised 3% of the total counts in the PIN band. 

Spectral fits used data from XIS 0 and 3, in the energy range 
0.75 - lOkeV and from the PIN over 15 - 50keV. Data in the 
range 1.75-1.9 keV were excluded from the XIS spectral analysis 
owing to uncertainties in calibration around the instrumental Si 

' http://www.astro.isas.jaxa.jp/suzaku/doc/suzakumemo/ 
suzakumemo-2007-01 .pdf 
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Table 1. Summary of XMM-Newton observations for Ark 564 


OBSID 

Start date 

Duration 

Instrument 

Net exposure 

Excision 




pn 

38.6 

8 

0670130201 

05/24/2011 

59.5 

RGSl 

55.6 

/ 




RGS2 

55.5 




pn 

34.3 

0 

0670130301 

05/30/2011 

55.9 

RGSl 

50.2 

/ 




RGS2 

50.1 




pn 

31.1 

0 

0670130401 

06/05/2011 

63.6 

RGSl 

43.0 

/ 




RGS2 

42.9 




pn 

40.6 

5 

0670130501 

06/11/2011 

67.3 

RGSl 

59.7 

/ 




RGS2 

59.6 




pn 

33.5 

0 

0670130601 

06/17/2011 

60.9 

RGSl 

48.4 

/ 




RGS2 

48.2 




pn 

29.1 

0 

0670130701 

06/25/2011 

64.4 

RGSl 

42.1 

/ 




RGS2 

41.9 




pn 

40.5 

4 

0670130801 

06/29/2011 

58.2 

RGSl 

58.1 

/ 




RGS2 

58.0 




pn 

38.7 

5 

0670130901 

07/01/2011 

55.9 

RGSl 

55.7 

/ 




RGS2 

55.7 


Notes. Observation ID, start date, duration of the observation (in kiloseconds), net exposure time after background flaring removal (in kiloseconds; 
the pn values include the loss of exposure due to the small window observing mode live time), and excision radius used to correct for pile-up (in 
arcseconds). 


K edge. XIS 1 was not used owing to the higher background level 
at high energies. The data were binned to have at least 20 counts 
per bin. In the spectral fitting, the PIN flux was scaled by a factor 
1.18, which is appropriate for the cross calibration of the XIS and 
PIN for this observation epochal A systematic error of 3% was 
added to the data. 


3. Spectral analysis 

We used Xspec v.12.7.0 for the spectral analysis. We used thex^ 
statistics as a goodness-of-fit measure, in a valid regime assured 
by the high number of counts in each energy channel for each 
analysed spectrum. All the mo dels include a Galac tic column 
density Nh = 5.43 x 10^° cm“^ (iKalberla et al.ll2005). modelle d 
with tbabs for the low-resolution pn data ( Wilms et al.ll200Cih . 
with tbnew for the high-resolution RGS data (Wilms et al. 2015, 
in preparation). Errors are given at the 90% confidence level 
(Ax^ - 2.706 for one parameter of interest). A cosmological 
redshift z = 0.02468 was adopted for the source. We first anal¬ 
ysed the RGS data to constrain the physical properties of any 
possible soft-band absorbing or emitting component (Sect. [ID; 
the results obtained for the soft band were then used as input 


ftp://legacy.gsfc.nasa.gOv/suzaku/doc/xrt/suzakumemo-2008-06.pdf 


in our modelling of the full-band EPIC pn and Suzaku XIS/PIN 
spectra (Sect. IjSl l. 

3.1. Grating analysis 

The RGS 1 and RGS2 stacked spectra were simultaneously fitted 
over the 0.4 - 1.9 keV band. Bad channels were ignored, and a 
systematic error of 3% was added in quadrature with the photon 
shot noise to account for the calibration uncertainties between 
the two RGS (see e.g. lKaastra et al.ir201 ll) . 

The spectral shape in this energy range is extremely steep, 
and a simple power-law emission model gives a very poor fit to 
the data (x^/v - 4700/1297). Adding a blackbody component 
significantly improves the fit statistic (Table ID. The black hole 
mass of Ark 564 is quite low (Mbh ~ 2.6 x 10^/V/^. lBotte et al] 
ItooI . so we might expect the tail of the thermal emission from 
the accretion disk to emerge in the soft X-ray band - but a simple 
blackbody emission component is obviously only a phenomeno¬ 
logical parametrization for the continuum shape. 

We replaced the power law -i- blackbody model with 
optxagnf, which represents the expected emissio n spectrum 
from the inner accretion flow around a black hole (iDone et al.l 
I 2 OI 2 I) . This model assumes that the accretion luminosity 
(parametrized in terms of the mass accretion rate Macc and the 
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Table 2. Fit to the RGS continuum 


Model 



Parameters 


Statistics 

bbody + po 

F 

0 ’7Q+0.02 

Nlkevl^O 

1 - 70 + 0.02 
^°- 0.02 

kTleW 

mil 

Nbh/lO-^' 

O+0.2 

■^•^-0.2 

xTy 

2596/1295 

optxagnf 

F 

2 48+0.13 

fpi 

0 60+° '® 

T/ LecIcI 

1 10+0.07 
^•^^- 0.02 

kT/eW T 

168:1^ 24^4 

Dorlrg x^ly 

80:®^ 2594/1293 


Notes. The power-law normalisation Nikev is listed in units of photons keV * cm ^ s * at 1 keV; the hlackbody normalisation Nt,b = L^g/D^^^, 
where L 39 is the luminosity in units of 10 '^^ erg s“* and Dio is the source distance in units of 10 kpc. 


black hole spin a) is produced in a disk + corona system. A geo¬ 
metrically thin and optically thick accretion disk extends from an 
outer radius up to the innermost stable circular orbit around 
the black hole and emits as a colour-temperature corrected black- 
body (i.e. it has an effective temperature higher than a standard 
SS disk at small radii within the region where electron scatter¬ 
ing starts to dominate the disk opacity: this effect is strongest 
for the hottest accretion disks, such as for low black hole mass 
and high mass accretion rate AGN) up to a radius rear- Inside 
this radius, the disk seed photons are Compton-upscattered in a 
two-phase hot medium: a fraction of them fpi is Comptonized on 
a very hot and optically thin medium (e.g. the “X-ray corona”), 
while the remaining (1 - fpi) is Comptonized on a much colder 
and optically thick medium (e.g. the “inner disk” or a “disk 
atmosphere” or a “disk skin”). We fixed the black hole mass 
Mbh — 2.6 X 10®Mo (iBotte et al.l l2004ll . the black hole spin 
a = CQ and the accretion disk outer radius rout — The free 

parameters are then the mass accretion rate Mace, the outer radius 
of the Compton-upscattering region rear, the temperature kT and 
optical depth r of its cold and optically thick portion, the frac¬ 
tion fpi of the energy that is upscattered in its hot and optically 
thin portion, and the observed photon index T. The optxagn 
model is also able to improve the fit statistics with respect to 
the power-law component alone and gives a statistically equiva¬ 
lent fit to the data as the power law -I- blackbody model (Ta¬ 
ble |2]l; however, given the two more free parameters required by 
optxagn and given the much faster performance of power law 
-I- blackbody, we proceeded to search for spectral residuals us¬ 
ing the latter as the baseline model for the continuum. First, the 
residuals were modelled using Gaussian emission or absorption 
lines (Sect. IrTTTT l. then a fit using Xstar tables was performed 
fSect. lTra . 

3.1.1. Fit with Gaussians 

A threshold of 99.9% in statistical significance was adopted to 
consider the addition of a Gaussian line to the model as required 
by the data, that is, > 16.3 for three degrees of freedom 
associated with each line: centroid energy, width, and inten¬ 
sity. The analysis was simultaneously performed on the RGSl 
and RGS2 dataset, and only features significant in both detec¬ 
tors were included in the model (except for energy ranges with 
missing chips, i.e. E - 0.90 - 1.18 keV for the RGSl, and 
E = 0.52 - 0.62 keV for the RGS2). 

We detected 16 absorption and 3 emission lines in the soft 
X-ray spectrum of Ark 564: they are shown as residuals against 

* The spin parameter a is degenerate with the mass accretion rate un¬ 
less one can constrain the latte r with UV flux measurements from the 
disk (see e.g. iDone et ^l2013l) . which we lack because of strong con¬ 
tamination of the host galaxy emission in this band; therefore we fixed 
for simplicity a = 0. 







Energy (keV) 

Fig. 1. RGS spectral residuals in the source rest frame. 


the continuum in Fig. [T] while their properties are listed in Ta- 
ble[3] For 14 out of 16 absorption lines we were able to associate 
a theoretical identification with H-like and He-like transitions of 
ions of C, N, O, Ne, and L-shell transitions of Fe: all of them 
are compatible with a zero-velocity shift. All the identified ab¬ 
sorption lines are narrower than the instrumental RGS resolution 
except for the features observed between ~ 7 00 - 800 eV , whic h 
correspond to the Fe I-XVI M-shell UTA (iBehar et al.l 1200 1|) . 
and the two high-ionisation transitions identified with Fe XVII 
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and Ne IX, which display a FWHM ~ 1000 - 1500 km s '. 
The three emission lines are identified with H-like transitions of 
N VI, N VII, and O VII. 

Two features remain unidentified at rest-frame energies of 
~ 552 eV and ~ 619 eV; both of them are resolved by the RGS, 
but do fall in regions that are devoid of known strong spectral 
transitions of cosmically abundant elements. We note that within 
the measurements errors, they might be identified with N VII and 
O VII Ka absorption by allowing a net blueshift of ~ 0.1c. 



Energy (keV) 




Fig. 2. Rest-frame plot of the three zones of warm absorbing gas 
detected in the RGS data. From top to bottom, the ionisation param¬ 
eter and column density are log(^/erg cm s“*) ~ -0.8, 1.3,2.4 and 
log(A;//cm-2) ~ 20.1, 20.3, 20.9. ” 


3.1.2. Fit with Xstar tables 

We generated a table using the latest version of Xstar 
V. 2.2 . Ibn3 and a finely spaced energy grid, using 10^ energy 
bins over the 0.1 - 20 keV energy range used for the compu¬ 
tation. The incident ionising continuum was a power law with 
a steep photon index T = 2.5, appropriate for Ark 564. So¬ 
lar abundances and a turbulent velocity of 100 km s ' were 
adopted. The gas column density was sampled over 10 points 
from Nh - 5 X 10'* cm“^ to Nh - 10^' cm“^, while the gas 
ionisation state was sampled over 20 points from log(^/erg cm 
s“') = -1.5 up to log(f/erg cm s“') = 3.5. 

The continuum emission was modelled with a power law 
plus a blackbody, and the model included the three emission 
lines detected in the previous section. We obtained a fair fit to 
the data = 1.33) using three zones of photo-ionised gas. We 


Table 4. Fit to the RGS spectra with Xstar 


Parameter 

F 

NlkeV 

kT 

Nbb 

logA^H.l 

log^i 

10gA^H,2 

log ^2 
log Ah, 3 
log ^3 


Value 

Tumr 


2.81_ 

1.81+' 


03)2 

+2 


+oh 


,+(?08 


136 
3.7 

20 . 10 _ 
-0.77+'^'^^ 
20 35 +^-^ 
134+0°'^' 
20 . 86 + 
2-36:o.o6 
1770/1331 


Notes. Units are keV^' cm“^ s“' at 1 keV for the power law normalisa¬ 
tion N\kev\ eV for the blackbody temperature kT ; cm“^ for the column 
densities NhT, erg cm s“' for the ionisation parameters f,. The black¬ 
body normalisation Nbb = where L 39 is the luminosity in units 

of 10 ^® erg s“' and Dio is the source distance in units of 10 kpc. 


found no net velocity shift for any of the three zones. The rela¬ 
tive contribution of each zone in terms of absorption components 
is shown in Fig.|2] The least ionised zone is responsible for most 
of the opacity for lowly ionised oxygen, that is, O IV, O V, and 
O VI; the intermediate zone reproduces the H-like and He-like 
transitions of C VI, N VI-VII, and O VII-VIII, as well as the 
least-ionised portion of the Fe M UTA; finally, the most-ionised 
zone contributes to filling the O VIII H-like absorption as well 
as the higher ionisation states of the Fe M UTA, the Fe L, and 
the H-like and He-like Ne lines. Fit parameters are reported in 
Tableg] 

We checked for partially covering gas along the line of sight 
by allowing any of the three absorbing zones to be only partially 
covering the source. While there is no statistical improvement 
when the zones are individually allowed to be partially covering, 
a - 30 for one extra degree of freedom is found when all 
the three zones share the same covering fraction C/ = 0.76+gQ2. 
While the gas ionisation state is not affected by introducing an 
unabsorbed component along the line of sight, the column den¬ 
sities measured in this scenario are of the order of 40% higher 
than those measured in the scenario Cf = 1. 


3.2. Broad-band analysis 

We show in Fig. |3] full-band (0.4-10 keV) pn light curves for 
each observation, binned to 128 s. The count rates for the four 
excised observations have been rescaled to account for the loss of 
counts due to discarding the central region of the PSF. Given the 
very low residual background count rate (typically lower than 
1% of the source one) after screening against the presence of 
strong flares, we show the source plus background light curves. 
The time variability of Ark 564 is striking on both short (kilosec- 
onds) and long (hours and weeks) time scales. We first present 
the spectral analysis of the time-averaged spectra (Sect. 13.2.1b 
and of the flux-selected spectra (Sect. 13.2.2b : we then present 
results of time-resolved spectral analysis fSect. lTXJb . 


3.2.1. Time-averaged spectra 

The eight pn spectra were stacked using mathpha. We gener¬ 
ated an average background spectrum by appropriately weight- 
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Table 3. Absorption and emission lines detected in the RGS stacked spectra of ARK 564 


Frest/eV 

crleV 

FWHM/km s' 

427 3+° "^ 

430.5+°] 

435 . 7 JJ 

0 9+°-^ 

< 0.2 
< 0.5 

1490+-^^° 

< 330 

< 810 

502.6^ 

< 1.2 

< 1690 

544.8^ 

552.1+°-® 

554.4+!^ 

< 0.6 

0 0 + 1.2 

< 0.4 

<780 

48601^540 

<510 

561.2+1 

< 0.5 

<630 

563.0^ 

< 0.5 

<630 

573.4^ 

< 1.0 

< 1230 

619.4® 

652.4+°-J 

2 0 +’-' 

z,.u_i 3 

< 0.5 

2280![2®o 

<530 

665.9+°-2 

745.8![°;® 

766.3^11 

812.9+: 

871.2® 

921.3/;°-' 

1019.7+^-° 

< 0.5 

74 0+'®-^ 
p 4 + 4.6 

< 2.2 

<530 
70100/;|4oo 
2210+|240 
< 1910 

1 1+°-'* 

+ ■ + - 9.5 

2 1+'-' 

^■ + - 1.3 

<4.4 

890+33° 

1610+*'^° 
iOiU JQQQ 

< 3050 


EW/eV 


0 34+'' “ 

-0^02 

-0 

-0 35“° ® 
0 31+° °8 

-0 89“^ " 
-1 43-0.18 
^■^■’+ 0.46 
-0 77-0 '0 
'+ 0.14 

-0J3 

-0 33“° '^ 
+fi4Q 

-0-48;: 

-0 60“°4o 
-069“^:?^ 

'^■'^^+042 

-13 

^-^-'^+ 3.20 

-0 72“°-^' 

_ 062 ^ 8 :i 

-1 29 “°'^ 

^ ■^^+0.09 
-1 37 "°'^' 


30.6 

83.3 

41.4 

25.5 
104.2 

47.3 

45.8 OV 1^-271 

55.1 OVII 

18.5 OVIli-2/7 

32.1 O VII ls- 2 / 7 (r) 

17.8 ? 

31.1 O VIII li-2/7 

70.6 O Villi-3/7 

170.5 Fe M UTA 

20.2 Fe M UTA 

17.4 Fe XVII 2/7 - 3<i 

75.8 Fe XVIII 2 / 7 -3t^ 

53.8 Ne IX Is- 2/7 

20.1 NeXls- 2/7 


430 (r), 426 (i), 420 (f) 
431 
436 
500 
545 
? 

555 

574 (r), 569 (i), 561 (f) 
564 
574 
? 

653 

666 

700 - 800 
700 - 800 
812 
873 
922 
1022 


ID 


Fth/eV 


N VI 

N VI Is-2/7 
C VI Is-3/7 
N VII Is-2/7 
OIV Is-2/7 


Notes. Col.(l): rest frame line centroid energy; Col.(2): observed frame line width; Col.(3): observed frame line width; Col.(4): observed frame 
line equivalent width; Col.(5): statistical improvement after adding the line, for three free parameters; Col.(6): line ID; Col.(7): theoretical line 
transition energy. All the errors are at the 90% confidence level. 


ing each individual background spectrum by the BACKSCAL 
ratio between each source and background extraction region, 
which varies among the different observations because of the 
different excision radii used to clean the spectra from pile-up 
effects. Average R MF and ARF were created using addnnf and 
addarf. Following lMiller et al.l (l2008l) . each individual spectrum 
was binned at the HWHM instrumental energy resolution. A sys¬ 
tematic error of 1% was added to the data below 2 keV to ac¬ 
count for uncertainties in the instrument calibration 0. Because 
of the different count losses due to excision of the PSF in differ¬ 
ent observations, we were unable to recover the intrinsic average 
flux of the source, hence we do not cite normalisation values for 
the spectral fit referred to the stacked pn spectrum. As a sanity 
check, we also fit simultaneously the eight spectra of each indi¬ 
vidual observation; in this case, we can recover the intrinsic flux 
of the source by correcting the count rate for the PSF count loss. 

The broad-band pn and XISh-PIN spectra are shown in Fig.|4] 
unfolded against a power law with T = 2. The spectra are ex¬ 
tremely steep and show some structure in the Fe K band. We 
first fitted the 3-10 keV spectra to a simple power-law con¬ 
tinuum emission. The fit statistics are^f^/v =484/512 for the 
XIS spectrum, = 190/113 for the stacked pn spectrum, and 
jv - 1084/905 for the simultaneous fit to the eight individual 
pn datasets. Spectral residuals to this model are plotted in Fig.|5] 
where the power law was fitted ignoring the Fe K region: even if 
a power law formally yields an excellent fit to the XIS data, sig¬ 
nificant residuals are evident. Adding a Gaussian emission line 
significantly improved the fit for all the datasets. Spectral pa¬ 
rameters for this model are reported in Table |5] The emission 
line centroid energy and width {Eune ~ 6.6 keV, crime ~ 400 eV) 
were found to remain constant within the errors in the different 

^ There are typically >>10^ counts per bin below 1 keV, ~ 10^ counts 
per bin in the Fe K hand, and ~ 500 counts per bin at the highest prohed 
energies. 



Energy (keV) 



Energy (keV) 

Fig. 4. Observed frame EFe XMM-Newton EPIC pn (top) and Suzaku 
XIS+PIN (bottom) spectra unfolded against a power law with F = 2. 

observations. The line equivalent width EW/,„(. was found to vary 
during the 2011 XMM-Newton observational campaign; it was 
higher when the source flux was lower. We return to this point in 
Sect. l3.2.^ when we analyse the flux-selected pn spectra. 

The single G aussian line was r eplaced by the reflection 
model xillver (iGarcla et al.l 1201 3ll . which simulates the re¬ 
flected spectrum of a power-law incident on a plane-parallel op¬ 
tically thick slab. We used solar abundances and tied the slope of 
the incident continuum to that of the primary power law. Good 
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Fig. 5. Rest frame spectral residuals for a fit to a power law in the 3-6 
and 7-10 keV band for the pn (top) and XIS (bottom) data. 


fit statistics (respectively;^^^/v =459/510,112/111, and 957/903 
for the XIS, the stacked pn, and the individual pn spectra) were 
obtained with a power law with F ~ 2.4 illuminating a highly 
ionised slab with log(^/erg cm s“') = 3.2 - 3.8 (90% c.l.). Al¬ 
lowing the normalisation of the reflected component to be a free 
parameter in the simultaneous fit to the eight individual pn spec¬ 
tra did not improve the fit. 

We then turned our attention to the full-band spectra, which 
were modelled including the three zones of photo-ionised gas 
and the three narrow emission lines detected in the RGS spec¬ 
tra (Sect. 13.1.21 1. There is a prominent excess of emission both 
at lower energies (down to £ ~ 0.4 keV) and at higher ener¬ 
gies (up to 30 keV with Suzaku). The broad-band spectral shape 
of the high S/N stacked pn spectrum can be roughly modelled 
by a combination of xiliver and either a power law plus a 
phenomenological blackbody (x^/v — 607/181), or optxagn, 
where the black hole spin was fixed to zero and the disk outer 
radius to lO^r^ (x^lv - 432/179). However, a simple combi¬ 
nation of optxagn plus an ionised reflector, as found for the 
fit to the 3-10 keV data, is unable to give a good fit to 
the broad-band data and leaves systematic residuals at £ > 5 
keV. A further cold reflection component (modeled again with 
xi liver) is able to well reproduce such high-energy residuals, 
but it strongly over-predicts the strength of soft emission lines at 
£ < 2 keV. The addition of a layer of ionised gas with a column 
\og{NH!cvcT^) ~ 23.5 in front of the cold reflector gave a good fit 
to the broad-band data by excluding such strong emission lines 
(see Table| 6 ]l. 

The same model (optxagn -H ionised reflection -H cold and 
absorbed reflection) was able to well reproduce the XISh-PIN 
spectra without significant variations of the model parameters 
with respect to the stacked pn ones. The observed variability for 
the eight individual pn spectra is driven by changes in the con¬ 
tinuum amplitude. We found significant improvement of the fit 
statistics {kx^ll^v = 151/7 and A^^/Av = 76/7, respectively) 
when allowing the photon index T and the fraction of energy 
dissipated in the hot corona fpi to vary between the observa¬ 
tions; these variations were very weak, however, with a max¬ 
imum AT ~ 0.1 and A/^/ ~ 0.05. The best-fit model spec¬ 
tral parameters are reported in Table | 6 ] while the Suzaku and 
stacked pn spectra are plotted in Fig.| 6 ]unfolded against the best- 
fit model, together with data-to-model ratios. 

A blurred reflection scenario We tested a blurred reflection 
scenario on the stacked pn spectrum. The intrinsic continuum 
was modelled with a power law plus a blackbody; a reflection 
component modelled with xillver was convolved with the rel- 




Energy (keV) 

Fig. 6. EPIC pn (top) and Suzaku XIS-l-PIN (bottom) spectra unfolded 
against the best-fit model described in Table|^ along with data-to-model 
ratios. The intrinsic continuum emission is modelled with optxagn 
(dotted-dashed black line); two reflectors are modelled with xillver, 
a highly ionised (dashed blue line) and a cold, absorbed one (solid ma¬ 
genta line). The total model is plotted in red. 


ativistic smearing model kdblur (iLaoil[19911) . The emissivity 
index q was fixed to 3 (a value typical of isotropic illumination 
of the reflecting slab) and the iron abundance was fixed to so¬ 
lar. A good fit ix^ Iv - 197/179) is obtained with a power law 
with r = 2.45 + 0.01 reflecting onto a ionised (log(^/erg cm 
s“') = 2.98/;QQg) slab seen edge-on (inclination i = 89 + 1 deg), 
that extends close to a rapidly rotating black hole (with a slab 
inner radius r,„ < 1.85rj at 90% c.l.); the unfolded spectrum is 
plotted in the top panel of Fig. |7] 

Replacing the power law -t blackbody model with 
optxagn gave a good fit to the data as well (y^/v = 167/177) 
with the following parameters: an intrinsic X-ray continuum 
with a slope F = 2.43 + 0.01; a mass accretion rate Mace - 
0 . 80 /QQ 2 -^£rf(/; a coronal radius Vcor - a fraction fpi - 

0.58/;qq 5 of accretion energy comptonized in the hot optically 
thin X-ray corona, the remaining (1 - fpi) Comptonized in 
a colder {kT - 152/2 eV), optically thick (t = 32 + 5) 
medium. The blurred reflection component has a high ionisa¬ 
tion state (log(^/erg cm s“') = 2.91 /qq 3) and extends down to 
r,„ < 1.59r^. The fraction of reflected or direct flux is 0.38 /qj 2, 
0.46/°°^, and 0.43/°“ in the 0.4 - 1, 1 - 4, and 4 - 7.5 keV 
band, respectively. The corresponding unfolded spectrum is plot¬ 
ted in the bottom panel of Fig.|7] A comparison between the top 
and bottom panels of Fig. |7] shows the equivalent modelling of 
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Table 5. Fit to the 3-10 keV band with a power law plus a Gaussian emission line 


Dataset 

F 

^/me/kcV 

Crtinel^V 

hine 

EW,,„,/eV 

T^/v / 3 - 10 /IO " ergcm ^ s ^ 

XIS 


^■‘^'^- 0.10 

901 +2J4 

^■^^-91 

"T’o+TTr" 

^■^-0.7 

138“! 

460/509 

1 

^■^'^-O.OS 

pn stacked 

2 4 Q+ 0.01 
^•^^- 0.02 

fz ^o+0.09 

460+“° 

^+’'+-100 

1 6+°+^ 

+ ■"- 0.4 

111+30 

96/110 

/ 

201 




< 2.3 

< 121 


l-63!°“ 

301 

401 




2 0+°° 
^■+'- 0.7 

2-5+H 

171+“ 


1 06+ ■' 

+ ■'+"-0.08 

501 

601 

9 4^+0.02 
^•^'^- 0.02 

^■'^'^- 0.10 

^10+160 

-120 

2 5^^:^ 

^■^^-0.9 

157!| 

166+“ 

925/895 

1 07+8:8? 

+ ■-’'-0.08 

1.24+°-°! 

701 

801 

901 




2 9+^'^ 

Am 

< 1.9 

TTI+fb! 
■’^ + -84 
990+70 

< 108 


n 77+D:fJ8 

12448 : 8 ? 

4 > 8:?8 
+ ■^-+-0.09 

pn low 
pn med 

9 4^+0.02 

6 57+0 09 
- 0.10 

= 460 

3 -lil 

245!“ 

142+“ 

585/553 

0.87!°:“ 

i- 3 i!:p 

pn high 




< 3.0 

< 81 


1 q<r+0.07 
-^-^-0.14 


Notes. Iiine in 10 ^ phot cm ^ s * - Xspec model: powerlaw + gauss. 

Table 6 . Fit to the full band with optxagn, an ionised reflector, and an absorbed cold reflector 


Dataset 

F 

T/ l^Edd 

rcorirg 

fpi 

kTleW 

T 

log Of 

log Nh 

^ahs 

log ^cold 

X^lv 

XIS-i-PIN 

2.51 +0.02 

1 19+0.04 

^•^■^-0.18 

32 + 2 

0 79+0.04 

^•'■^-0.07 

140!! 

61+5 

3 46+u.^4 

^-0.07 

99 1 +0.6 

-0.3 

< 2.9 

< 1.6 

1090/1109 

pn 

9 cq+0.02 
^*^^-0.01 

^•^^-0.13 

> 35 

0 75+0-02 
^•/-’_0.10 

148“ 

41 + 15 

‘+ + _7 

^■^^- 0.43 

99 47 + 0.10 
-0.19 

< 1.6 

0 . 86 “;“ 

196/175 


201 

301 

401 

501 

601 

701 

801 

901 

pn low 
pn med 
pn high 


2.63+° 

2.63 

2.57 
2.59 

2.58 

2.59 
2.54 
2.63_„ 


02 


02 


2 66 +° °^ 


1 46+°-^ 
+ ■^+’- 0,2 

1 . 02 +°“ 
1.03+ 
1.16+ 
1.03+ 
0 . 66 + 
0.99+ 
1.34+ 


0.65+ 

1.03+ 

1.60+ 


58!“ 


0.76!° 

0.75!° 

0.76!° 

0.75!° 

0.76+° 


02 


;8:85 

48:81 

-0.02 

0.71!| 

0 75+0-02 
^ ''^-0.03 

+0.05 


48:85 59+13 

48:8t 

n 89+0.03 

^■O-^_0.02 

-0.07 



150!^ 38!f 


144“ 


> 51 


3.5 


+ 1.2 
1.1 


23 


a Q1 +0.61 
*^•^^- 0.06 


2 10 +°“ 
^■ + '+- 0.16 


99 7J_+0.09 

/H_o 15 


< 0.52 1491/1457 


1.97!°“ <0.52 


1008/907 


Notes. Xspec model: optxagnf + xillver;,o, + absorber*xillvercoH. Units are cm ^ for the column density N//; erg cm s ' for the ionisa¬ 
tion parameters 


optxagn and a power law -i- blackbody where the peak tem¬ 
perature of the latter is kT ~ 130 eV. 

The application of the blurred reflection scenario to the 
broad-band Suzaku data, with an emissivity index fixed to ^ = 3 
and iron abundances fixed to solar, gave a good fit (y^/v = 
1084/1109) with spectral parameters compatible within the er¬ 
rors with those derived for the stacked pn spectrum, but for the 
inner radius of the reflecting slab, which is found to be much 
higher; = 27!2“rg. The XISh-PIN spectrum unfolded against 
this model is plotted in Fig. 0 despite the formally excellent fit, 
significant residuals in the PIN band are present, with the source 
flux being underestimated by the blurred reflection model. 

A partial covering scenario A partial covering scenario was 
tested by absorbing the intrinsic continuu m (modelled with 
optxagn) with zxipcf (iReeves et alJl2008l) . Both the stacked 


pn spectrum and the XIS-i-PIN spectra are well modelled by 
these two components (y^/v = 198/178 and 1095/1110, re¬ 
spectively). The absorber column density, ionisation parame¬ 
ter, and covering fraction were log/An/cm^^) = 23.62!g°2, 
log(f/erg cm s“') = and C/ = 0.43!°Qg for the pn, and 

log/An/cm^^) = 23.75!°'^°, logf/(erg cm s“^) = 1 . 9 !° 4 , 

Cf - 0.29!^ °^ for the XIS-hPIN data. The absorbed component 
helps reproducing the spectral shape at high energies, basically 
replacing the absorbed reflector plotted as the magenta line in 
Fig.| 6 ] The intrinsic continuum model parameters were similar 
for the two datasets and comparable with the parameters derived 
in the previous section and reported in Table | 6 ] except for the 
mass accretion rate that in this scenario is estimated to be higher, 
Macc/MEdd = and Macc/MEdd = 1.51!°“ for the pn 

and XIS-hPIN data. 


Article number, page 8 oflH 

















M. Giustini et al.: X-ray reverberation within one light hour of the central SMBH in Ark 564 


Table 7. Flux-selected spectra 


Flux state 

OBSID 

Exposure 

Count-rate 

Stacked set 

Counts(4-10) 


301 

11.4 

26.26 + 0.05 



Low 

401 

601 

10.3 

7.6 

25.97 + 0.05 
26.71+0.06 

301+401+601+701 

25600 


701 

24.8 

20.41+0.03 




301 

9.5 

37.05 + 0.06 




401 

9.4 

37.7 + 0.06 

301+401+601 

18800 

Medium 

601 

8.2 

37.21+0.07 




501 

19.8 

38.24 + 0.04 

501 

9200 


801 

23.2 

37.88 + 0.03 

801 

11900 

High 

201 

901 

12.4 

9.7 

57.16 + 0.05 
61.41+0.05 

201+901 

10300 


Notes. Notes: flux state name, OBSID, exposure (including the live time), count-rate in the 0.4-10 keV band (corrected for the psf excision), 
stacked spectra, net counts in the 4-10 keV band. 



_1_1_I_ . . I _1_I_1_L. 
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Energy (keV) 
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Fig. 7. Stacked pn spectrum unfolded against two different represen¬ 
tations of a blurred reflection scenario. Top panel: continuum modelled 
with power law (black dotted-dashed line) + blackbody (blue dot¬ 
ted line), and blurred reflection modelled with kdblur*xillver (blue 
dashed line); bottom panel: same as above, but modelling the continuum 
with optxagn (black dotted-dashed line). 


3.2.2. Flux-selected spectra 

The EPIC-pn source spectra were stacked according to three dif¬ 
ferent flux levels selected after inspecting Fig. |3] a low flux state 
with a count rate < 27 counts s“*, a medium flux state with 



Energy (keV) 

Fig. 8. Suzaku XIS-l-PIN spectra unfolded against a blurred reflection 
model, colour-coded as the bottom panel of Fig. [T] 



Energy (keV) 


Fig. 9. 0.4 -10 keV flux-selected spectra unfolded against a power law 
with F = 2: low (black), medium (red), and high (blue) flux states. 

30 counts s ' < count rate < 40 counts s ', and a high flux 
state with a count rate > 45 counts s“'. Flux-selected spectra 
were then stacked according to the shape of the source extrac¬ 
tion region and were fitted if they had > 5000 net counts in the 
4-10 keV band. The relative contribution of the different ob¬ 
servations to each flux-selected spectrum is reported in Table [T] 
while the full-band unfolded spectra are plotted in Fig.|2 despite 
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flux changes of over a factor of two, the spectral shape does not 
vary dramatically. 

The 3-10 keV band of the three flux-selected spectra was 
fitted simultaneously to a power-law model; the fit statistics was 
quite poor (x^ jv - 708/557), and there were residuals in the Fe 
K band. As with the singular spectra, the addition of a Gaussian 
emission line (where the width cr was fixed to 460 eV, i.e. the 
best-fit value found for the average stacked pn spectrum) signif¬ 
icantly improved the fit (A;^f^/Av = 123/4). The line centroid 
energy is consistent with the one measured in the average spec¬ 
tra, but its equivalent width significantly decreases from lower to 
higher states (see Table|5]for values). This suggests the presence 
of a constant (or varying slower than the continuum) component 
that emerges when the intrinsic continuum is in a low state. 

We applied the best-fit model found in Sect. l3.2.TI (optxaan 
+ ionised xillver + absorbed cold xillver) to the flux- 
selected spectra, which yielded a good fit to the data (x^/v - 
1008/907, Table |6]) with the two reflection components remain¬ 
ing constant, and the observed variability being driven by varia¬ 
tions in intrinsic continuum amplitude. 

3.2.3. Time-resolved spectra 

We aim at extracting the spectral signature of the ’hard delayed 
ex cess’ that follows flares of emission in Ark 564, as identified 
bv iLegg et al.l (l2012h . To this end, we extracted the source spec¬ 
tra in time intervals centred on the position of the ’soft flares’B 
with a width of +384 s. The associated delayed signal was ex¬ 
tracted in a time interval 1088 - 2496 s after each peak centre, 
that is, the time interval when t he hard delayed e xcess appears 
to be the strongest (see Fig. 6 of iLegg et al.ll2012l) . All the ’soft 
flares’ spectra were then stacked together, as were all the ’de¬ 
layed’ spectra. For the spectral analysis we used the same scaled 
background as we used to analyse the average stacked pn spec¬ 
trum. Figure [10] shows the unfolded broad-band spectra for the 
flaring (black) and their delayed (red) time intervals. The vari- 
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Fig. 10. Top panel: soft flares (black) and their associated delayed sig¬ 
nal (red) spectra unfolded against a power law with F = 2. Bottom 
panel: data-to-model ratio to a power law fitted on the 2-10 keV band. 


^ Soft flares are defined bv iLegg et all ( 120121) as those only detected in 
the 0.4 - 1 keV band. 


ation in spectral shape is clear: while at energies E < 2.5 keV 
the spectral shape is just shifted down in intensity (as found for 
the flux-selected spectra), at £ > 3 keV there is an evident hard¬ 
ening of the delayed spectrum. This can be seen in the bottom 
panel of the same figure where high-energy residuals to a fit to a 
power law with F ~ 2.4 are plotted. 

We modelled the two spectra exploiting t he con straints de¬ 
rived from the timing analysis bv iLegg et alJ (1201 2h in terms of 
relative fractions of direct (flaring) and scattered (reprocessed) 
light. In particular, as the low-frequency lags are always positive 
with respect to the 0.4 - 1 keV band, the fraction of scattered 
light must increase from the softest to the hardest band: we used 
this constraint to construct a model able to reproduce the flaring 
and delayed spectra. 




Energy (keV) 


Fig. 11. Top panel: soft flares (black) and their associated delayed 
signal (black) spectra unfolded against a reflection-dominated model: 
optxagn (dotted-dashed lines), a heavily absorbed cold reflector (dot¬ 
ted lines), and a slightly absorbed hot reflector (dashed lines). Bottom 
panel: same as above, but unfolded against a Compton-upscattering- 
dominated model: optxagn (dotted-dashed lines), a hot reflector (blue 
dashed line), and a hot Compton-upscattering model (dotted lines). 

First, we tested the best-fit model of Sect. l3.2.TI (i.e. optxagn 
plus an ionised reflector plus a cold, absorbed reflector, both 
modelled with xillver). Within this model, optxagn repre¬ 
sents the direct light, while the two reflected components rep¬ 
resent the scattered or reprocessed light. The normalisation of 
the continuum was allowed to vary between the two data sets 
and is ~ 10% higher in the flares than in the delayed spec¬ 
trum. To match the constraints on the scattered or direct light 
fractions in the different energy bands, another layer of ionised 
gas was introduced in front of the ionised reflector. This model 
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can reproduce the data very well (x^ jv - 350/355) with con¬ 
tinuum emission model parameters compatible with those mea¬ 
sured for the average spectra (Sect. [JS/Tl Table | 6 ]l, but a higher 
flux in the two reflected components for the delayed signal; the 
0.4-10 keV observed flux is l-h/gg x 10“'' erg s“' cm“^ and 
X 10“'^ erg s“' cm“^ for the ionised and cold reflector in 
the delayed spectrum, respectively, as opposed to 1 .3/;Qg x 10 “" 
erg s“' cm“^ and ^ x 10 “'^ erg s“' cm“^ in the flaring spec¬ 
trum. The ionised reflector is absorbed by a layer of gas with 
log(AH/cm“^) = 21 . 6 /;q j and log(^/erg cm s“') = 0 . 21 /;° [^. 
The corresponding model and spectral residuals are reported in 
the top panel of Fig. [TT] 

Alternatively, reprocessing may occur onto a hot medium 
such as to be able to Compton-up-scatter the input (soft flares) 
photons. We tested this scenario by replacing the co ld absorbed 
reflector with the comptt model (lTitarchuHll994h . where we 
assumed a disk-like geometry and a temperature for the seed 
photons equal to the temperature of the warm, optically thick 
medium responsible for the Compton-upscattering of disk pho¬ 
tons inside rear as derived from optxagn, that is, ~ 

150 eV. Within this model, optxagn again represents the di¬ 
rect light, while the scattered or reprocessed light is repre¬ 
sented by comptt. This model can also well reproduce the data 
(X^/v - 380/360) by allowing again the optxagn normalisa¬ 
tion to be higher in the flare spectrum than in the spectrum with 
the delayed signal, and the Compton-upscattering component to 
emerge in the delayed spectrum, with a 0.4-10keV observed flux 
of 5 . 6 /;q 3 X 10 “'^ ergcm“^ s“' and 6 . 6 /q 2 X 10 “'^ ergcm“^ s“' for 
the flares and delayed spectrum, respectively. The temperature of 
the Compton-upscattering material is found to be kT^ - 4.3/q ^ 
keV, its optical depth t*' = h/j- The flares and delayed spec¬ 
tra unfolded against this model, as well as spectral residuals, are 
plotted in the bottom panel of Fig. [TT] 




TIME s Energy (keV) 


Fig. 12. Left panel: a segment of the 0.4-10 keV light curve of the pn 
observation 901 where a strong and sudden increase and decrease in 
flux is evident. Right panel: four 0.2-10 keV EPIC-pn spectral slices 
extracted at the beginning (during the first ks) of the high-flux level 
(black), during the remaining of the high-flux level (red), during the first 
ks of the decrease in flux (green), and during the remaining decreasing 
flux level (blue): the time intervals are marked in the top panel with the 
same colour-code. Error bars have been omitted for clarity. 

Detailed flare evolution We have shown that the rapid source 
variability is dominated by intrinsic continuum amplitude vari¬ 
ations (Sect. 13.2.11 1 and that a hard spectral component appears 
delayed with respect to the peaks of the emission (Sect. 13.2.31 1. 
A plausible interpretation is that the continuum light reverber¬ 
ates on some reprocessing material located ~ 1000 light-seconds 
away from the continuum source. If this is the case, we would ex¬ 


pect to detect the reaction of this reprocessing gas following the 
longer-term source variability. 

To test this expectation, we selected a time interval during 
the XMM-Newton observation 901 where there is a strong long¬ 
term flaring emission with large amplitude, namely the 7000- 
17500 s segment. The total count rate undergoes amplitude os¬ 
cillations of about a factor of two on a time scale of a few kilo- 
seconds. In particular, there is a strong rise in flux just after the 
first 5 kiloseconds of the observation, when the flux increases by 
50% in about 1500 s. The flux remains on a high level for about 
6 ks, then decreases by a factor > 2 in about 5000 s (Fig.[T 2 | left 
panel). The corresponding spectra are plotted in the right panel 
of Fig. dll omitting the error bars for clarity; spectral variability 
in response to the flux changes is evident. We can parametrise 
the spectral hardness by means of the ratio R of the observed 
count rate in a hard (4.5 - 7 keV) and in a soft (0.4 - 1 keV) 
energy band. During the first kilosecond of the flare (black spec¬ 
trum) the source is soft, with R - 0.012; after one kilo second 
has elapsed, the hard band has had time to react to the increased 
soft flux (red spectrum) and is found to be in a medium hard¬ 
ness state, R - 0.017. Then the flux starts to drop; during the 
first kilosecond of decreasing flux (green spectrum) the hardness 
is highest (R - 0.019), given that the hard band still did not 
have time to adjust to the soft-band flux changes; then, after one 
kilosecond has elapsed, the source is found again in a medium 
hardness state (blue spectrum, R - 0.017). 


4. Results 

We present the results of both time-averaged and time-resolved 
X-ray spectroscopy of eight XMM-Newton and one Suzaku ob¬ 
servations of the NLS1 Ark 564. 

- The Ark 564 X-ray continuum appears extremely steep in the 
whole energy band analysed in this work (0.4 - 30 keV, see 
Fig. lUi. A fit to the RGS data in the range 0.4 - 1.9 keV 
gives an intrinsic slope F = 2.5 - 2.8 depending on the 
model adopted, optxagn or (powerlaw -i- blackbody) 
(Table|2]). Modelling the broad-band data with a more phys¬ 
ical model (optxagn -i- two reflectors) reveals a hard 
spectral component emerging above E > A keV; the in¬ 
trinsic continuum slope of Ark 564 thus ranges between 
F = 2.5 - 2.7 (Tableig. 

- High-resolution X-ray spectroscopy was performed in the 
soft band {E < 1.8 keV) with the RGS instrument (Sect. 
EB, and revealed the presence of ionised gas along the line 
of sight in the form of moderately ionised C, N, O, and Fe 
(Fig.djand TableO. 

- Modelling the absorption features with the 2D photoionisa¬ 
tion code Xstar requires three layers of absorbing gas, with 

ionisation parameter log(^/erg cm s“')-0.8, 1.3, and 2.4 

and no velocity shift with respect to the systemic redshift of 
the source; the total measured absorbing column ranges from 
~ 10^' cm"^ (fully covering absorber) to ~ 1.4 x 10^' cm“^ 
(allowing the absorber to be partially covering the continuum 
source), see Fig.|2|and Table |4| 

- A weak feature (EW~ 150 eV, see Fig.|5]l is detected in emis¬ 
sion in the Fe K band, and it can be modelled either with a 
single broad Gaussian (cr ~ 300 - 600 eV) or with three nar¬ 
row Gaussian emission lines with energies fixed to 6.4, 6.7, 
and 6.97 keV. The gas responsible for creating such features 
has a high ionisation state; replacing the Gaussians with a 
physically self-consistent reflection model such as xillver 
requires log(^/erg cm s“') > 3.5. 
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- Overall, the broad band spectrum of Ark 564 can be well 
reproduced by the superposition over the intrinsic continuum 
(modelled with optxagn) of two reflected components, both 
modelled with xiliver; one hot (log(^/erg cm s“') ~ 3-4), 
and one cold (log(^/erg cm s“^) < 1.2) and absorbed by a 
column density log(A///cm“^) ~ 23 - 24 (see Fig. |6] and 
TableHl). 

- A blurred reflection scenario was also tested on the time- 
averaged spectra by convolving the reflected component 
with a relativistic convolution model (i.e. kdblur). Both the 
EPIC-pn and the Suzaku spectra can be reproduced statisti¬ 
cally well in this scenario (Sect. 13.2.11 Fig. |7]), albeit giving 
discrepant results in terms of the inner radius of the reflecting 
slab, which would be r,„ < 1 .SSr^ in the case of the EPIC-pn, 
and rin = 27;|;2o^^.? the case of the Suzaku data. This model 
underestimates the Suzaku PIN flux (see Fig.[^. 

- A partial covering scenario is also able to well reproduce 
the EPIC-pn and Suzaku spectra, assuming that the intrinsic 
continuum is partially covered (C/ ~ 0.4 - 0.3) by a column 
of log(A///cm“^) ~ 23.6 - 23.7 of moderately ionised gas 
(log(^/erg cm s“t) ~ 1 - 2), for the EPIC-pn and the Suzaku 
XIS-hPIN data, respectively. 

- Given the highly variable X-ray flux of Ark 564 (see Fig.[3]l, 
spectral analysis was performed on three flux-selected spec¬ 
tra (Sect. [2221). Despite variation in flux of a factor > 2, the 
source spectral shape is does not vary dramatically (Fig. 0. 
Applying the best-fit model for the average spectra (optxagn 
-H hot xillver -t- cold absorbed xillver) to the flux- 
selected spectra yields a good fit to the data, with the ob¬ 
served variability found to be driven by intrinsic continuum 
amplitude variations. 

- Time-resolved spectral analysis was performed on the eight 
XMM-Newto n datasets b y stac king all the time intervals 
identified by iLegg et ^ (120121) as ’soft flares’ as well as 
all the subsequent ’delayed’ emission (Sect. 13.2.31) : the ’de¬ 
layed’ excess of emission that follows the soft flares is spec¬ 
troscopically identified as a spectral hardening at £ > 4 keV 
(Fig.®. 

- The best-fit model for the time-averaged spectra was applied 
to the ’soft flares’ and the ’delayed’ spectra by using the addi¬ 
tional constraints of the relative fractions of direct light (the 
’soft flaring’ signal) and scattered or reprocessed light (the 
’delayed’ signal) a t different energies , as measured from the 
timing analysis by iLegg et alJ (1201 2h to ’tune’ the different 
spectral components. 

- In particular, from the notion that the fraction of scattered or 
direct light must be decreasing from high to low energies, 
we added to this model another layer of ionised gas with 
log(AH/cm“^) ~ 21.6 and log(^/erg cm s“^) ~ 0.2 in front 
of the ionised reflector. The model is thus optxagn -i- hot 
absorbed xillver -i- cold absorbed xillver. Within 
this model the observed spectral variability is driven by a 
decrease in continuum amplitude and an increase of the re¬ 
flected flux in the ’delayed’ spectrum compared to the ’flar¬ 
ing’ spectrum. Given the time scales involved, a plausible 
physical scenario is to assume that the soft photons emitted 
during the continuum flare interact with the atoms of a gas lo¬ 
cated ~ lOOr^ from the continuum source and are scattered, 
which produces the hard ’delayed’ emission. 

- The gas responsible for scattering the soft flare emission 
could also be so hot as to be extremely ionised and thus 
be able to Compton up-scatter the flaring soft photons. We 
tested this scenario by replacing the cold absorbed reflector 
with the Comptonization model comptt and by removing 


the layer of absorbing gas in front of the ionised reflector 
(given the spectral shape of comptt, it is not necessary in 
this scenario). Again, we tuned the model following the con¬ 
straints of the relative fractions of direct or scattered light in 
different energy bands. We found that this model is also able 
to reproduce the data with respect to the scattered or direct 
light fraction constraints in all energy bands,if the gas has a 
temperature kT ~ A keV and a moderate optical depth t ~ 6 
(Fig. dB. 

- A qualitative analysis of the detailed evolution of a single 
strong flare was presented. We showed that the source spec¬ 
tral shape hardens ~ 1 ks after a rise in flux, and that con¬ 
versely, after a drop in flux it remains hard for ~ 1 ks before 
returning to its ’usual’ spectral shape (Fig.fTB. This is hard 
to reconcile with relativistically blurred reflection-dominated 
scenarios because we would expect the soft and hard band 
to ’communicate’ on much shorter time scales (correspond¬ 
ing to a few light travel times instead of to ~ a hundred), 
and strongly points, again, to reprocessing of the continuum 
photons by ionised gas located at lOs-lOOs from the cen¬ 
tral SMBH. 

5. Discussion and conciusions 

Ark 564 has a very steep X-ray continuum: a simple power-law 
model fit to the Suzaku and EPIC-pn data above 3 keV yields a 
slope of r ~ 2.4 (Table |5](. However, the intrinsic continuum is 
even steeper, being observationally diluted by a reflected or re¬ 
processed component that emerges in the hard band above £ > 4 
keV both in the time-averaged and in the time-resolved spec¬ 
tra. In a standard two-phase model for the X-ray emission of 
AGN (a cold optically thick accre t ion di sk -i- a hot optically thin 
corona. [Haardt & Maraschilll99lL [r993h . such a steep intrinsic 
continuum would suggest a high optical depth (0.5 < t < 1) 
and/or a low tempe rature (0 < 50 keV) for the hot X-ray corona 
(iHaardt et aklflo^ . 

To first order, the broad band spectra of Ark 564 are 
well modelled either by (power law -i- blackbody) or by 
optxagn and do not vary dramatically in shape despite flux vari¬ 
ations of the order of a factor of ~ 2 (Fig. 0. There are clear 
signs of reprocessing gas in the inner regions of Ark 564, both 
in the soft and in the hard X-ray band. 

In the soft band, a low column density of moderately ionised 
gas was found in the form of three different ionisation zones; the 
lowest ionisation state component of this absorbing gas might 
be the same as was o bserved in the UV band and studied by 
ICrenshaw et al.l (l2002h . who measured a modest velocity shift 
vuv ~ 200 km s“' and a column density Nh,uv ~ 1-5 x 10^^ 
cm“^. Three zones of ionised gas, with a comparable total col¬ 
umn density as measured in this work, as well as a modest 
(an upper limit of Vgut < HO km s~ *) outflowing veloc ity for 
the warm absorber was also found bv ISmith et al.l (l2008h in the 
analysis of five stacked XMM-Newton observations of Ark 564 
performed between 2000 and 2001. A Chandra/HETG observa¬ 
tion performed in 2000 also revealed multiple phases of ionised 
gas with a modest colu mn density and zero net velocity shift 
(iMatsumoto et al.ll2004h : new observations performed with the 
same instrument in 2008 res olved a velocity of a bout 100 km s“' 
for the multiphase absorber (iGupta et al.l2()T^ . Using the same 
dataset, iGupta etall (1201 3bl) claimed the discovery of a mildly 
relativistic {Vom ~ 0.1c) outflowing absorber by identifying three 
absorption features due to OVI and OVII. We note that the RGS 
data show two non-identifled absorption features (Fig.[T]), which, 
if identified with NVII and OVII Ka transitions, would suggest a 


Article number, page 12 ofll4l 




































M. Giustini et al.: X-ray reverberation within one light hour of the central SMBH in Ark 564 


similar outflow velocity Vout ~ 0.1c. However, the identification 
of a wind phase based solely on two absorption troughs where 
most reasonable photo-ionisation models would predict a wealth 
of absorption lines due to many species is challenging, and we 
do not discuss these features further. 

In the hard X-ray band, a weak feature (EW~ 150 eV) is 
detected in emission in the Fe K band, and it can be modelled 
either with a single broad Gaussian (cr ~ 300 - 600 eV) or with 
three narrow Gaussian emission lines with energies fixed to 6.4, 
6.7, and 6.97 keV. In any case, ionised iron in the inner regions 
of Ark 564 is required by the data. Given the width of the single 
Gaussian, simple Keplerian motion of a distribution of clouds 
around the SMBH would suggest a distance of the order of 100- 
1000s rg. The same range of distance would be suggested by 
the lack of reaction of the iron emission complex to the large 
amplitude variations of the ionising flux (Sect. lTZ^ . 

However, to reproduce the spectral shape at the highest en¬ 
ergies probed in this study, a hard-dominated component is 
required. This can be modelled either with a second (colder) 
reflector, a partially covered component, or with a relativisti- 
cally blurred reflector. In the cold reflection scenario, a column 
Nh > 10^^ cm“^ must absorb the reflected emission to not over¬ 
predict the strength of soft X-ray emission lines, which results 
in a quite contrived geometry. In the partial covering scenario, 
during the four years elapsed between the Suzaku and XMM- 
Newton observations the absorber covering factor would have 
significantly decreased while its ionisation state would have in¬ 
creased; in addition, the mass accretion rate in this scenario 
would be super-Eddington, which is unlikely. In the blurred re¬ 
flection scenario, the XMM-Newton data require an accretion 
disk extending down to < 2rg of the black hole, which is incon¬ 
sistent with the Suzaku data, which require r,„ ~ SOr^. Noting 
also how the inclination angle of ~ 90 degrees would be hard 
to reconcile with t he fact that Ark 564 is seen almost face-on 
dSmith et al.ll2008h . and the severe underestimation of the PIN 
flux in this sce nario (Fig.l^, we d o not consider this to be physi¬ 
cally plausible. iKara et al.l(l2013h reported the discovery of high- 
frequency Fe K lags of about 100 s in the first half of the 2011 
XMM-Newton observational campaign on Ark 564 and ascribed 
them to relativistic reflection close to the central SMBH. It is 
unclear, however, how to reconcile this scenario with the lack of 
any lag during the second half of the observational campaign. 
Relativistic effects might be at work in the inner regions of Ark 
564, but probably not that extreme: given the black hole mass 
estimate for Ark 564, a time l ag of 100 seconds corresponds 
to about ten gravitational radii. lEmmanoulopoulos et'afl (12014h 
inferred in the framework of a relativistic reflection lamp-post 
model a low spin of about 0.05 for the central SMBH of Ark 
564, which again suggests that relativistic effects, if present, are 
not dominant in the very inner regions of Ark 564. 

Alternatively, reprocessing of the soft photon flares could oc¬ 
cur in a medium so hot as to be able to Compton-upscatter them, 
giving rise to the hard excess component. In this case, we would 
have three different regions where Compton-upscattering occurs: 
a fraction (about 75%) of the UV seed photons from the accre¬ 
tion disk within ~ 5Qrg from the central supermassive black hole 
would be Comptonized in the very hot and optically thin X-ray 
corona (t << l,kT ~ 100 keV), giving rise to the primary X- 
ray power-law emission; the remaining ~ 25% of the UV seed 
photons would be Comptonized in a warm and optically thick 
medium (t ~ 40, kT ~ 150 eV) that contributes to the emission 
below ~ 2 keV; and all of this Comptonized emission would be 
then further Compton-upscattered in a hot and optically thick 
medium (t ~ 6,kT ~ 4 keV) located 10s to 100s of gravitational 


radii from the central supermassive black hole, following flares 
of emission. It would then follow that Comptonization plays a 
fundamental role in explaining the physics of the inner regions 
of AGN: in recent years, evidence is growing for it to explain 
the b road-band spectrum of several AGN, see, fo r instance, Mrk 
509 (iMehdipour et al.ll201 UlPetrucci et al.ll2013h. PG 1244-H026 


(|.lin et al.ll2013h. IH 0419-577 dPi Gesu et al.ll2014l). Ark 120 


dMatt et ai.l|20ri . and NGC 5548 (iMehdipour et al.lUOllh . and 
references therein. Future high-quality spectral and timing ob¬ 
servations on extended energy ranges with satellites such as 
NuSTAR (iHarrison et al.ll201^ will be able to distinguish be¬ 
tween the different physical scenarios and test the hypothesis that 
Comptonization is the main mechanism at work. 
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Fig. 3. Light curves for the eight XMM-Newton observations extracted in the 0.4 - 10 keV band and binned to 128 s. 
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